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We investigate the thermodynamic properties including equation of state, the trace anomaly, the 
sound velocity and the specific heat, as well as transport properties like bulk viscosity in the Z(2) 
and 0(4) models in the Hartree approximation of Cornwall- Jackiw-Tomboulis (CJT) formalism. 
We study these properties in different cases, e.g. first order phase transition, second order phase 
transition, crossover and the case without phase transition, and discuss the correlation between the 
bulk viscosity and the thermodynamic properties of the system. We find that the bulk viscosity over 
' entropy density ratio exhibits an upward cusp at the second order phase transition, and a sharp 

peak at the fst order phase transition. However, this peak becomes smooth or disappears in the 
case of crossover. This indicates that at RHIC, where there is no real phase transition and the 
system experiences a crossover, the bulk viscosity over entropy density might be small, and it will 
not affect too much on hadronization. We also suggest that the bulk viscosity over entropy density 
ratio is a better quantity than the shear viscosity over entropy density ratio to locate the critical 
endpoint. 
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I. INTRODUCTION 



: 

Studying Quantum chromodynamics (QCD) phase transition and properties of hot/dense quark matter at high 
temperature and baryon density has been the main target of heavy ion collision experiments at the Relativistic Heavy 
y—( . Ion collider (RHIC), the forthcoming Large Hadron Collider (LHC) and FAIR at GSI. 

^ ' At small baryon chemical potential /i, for QCD with two massless quarks, the spontaneously broken chiral symmetry 
is restored at finite temperature, and it is shown from lattice QCD [l[ and effective QCD models Q that this phase 
transition is of second order and belongs to the universality class of 0(4) spin model in three dimensions [|[. For 
| real QCD with two quarks of small mass, the second order phase transition becomes a smooth crossover at finite 
. ■ temperature. At finite baryon chemical potential, there are still no reliable results from lattice QCD due to the severe 
fermion sign problem. However QCD effective models Q suggest that the chiral phase transition at finite [i is of first 
order. It is expected that there exists a critical end point (CEP) in the T — fj, QCD phase diagram. The CEP is 
defined as the end point of the first order phase transition, and belongs to the Z(2) Ising universality class The 
precise location of the CEP is still unknown. In the future plan, RHIC is going to lower the energy and trying to 
locate the CEP. The signature of CEP has been suggested in Refs. Q. Recently, the authors of Ref. [|| suggested 
using the shear viscosity over entropy density ratio rj/s to locate the CEP . 

The ratio of shear viscosity over entropy density r]/s has attracted a lot of interests. It was expected that deconfincd 
quark matter formed at high temperature should behave like a gas of weakly interacting quark-gluon plasma (wQGP). 
It is now believed that the system created at RHIC is a strongly coupled quark-gluon plasma (sQGP) and behaves like 
a nearly "perfect" fluid @, Q- O ne crucial quantity is the shear viscosity over entropy density rj/s, which is required 
to be very small and close to the the lower bound r\js = 1/Aix Q to fit the elliptic flow at RHIC from hydrodynamic 
simulation [To| . Lattice QCD calculation confirmed that r\J s for the purely gluonic plasma is rather small and in the 
range of 0.1 — 0.2 [ll[. The perturbative QCD calculation gives a large shear viscosity in the wQGP with r\j s ~ 0.8 
for a s — 0.3 [TH]. Recent calculations using Boltzmann approach of multiparton scatterings (BAMPS) show that the 
small shear viscosity over entropy density can also be obtained by considering perturbative QCD inelastic scattering 
99 ->■ 999, sec Ref. [Hj]. 

In fluid dynamics, there is another important transport coefficient, the bulk viscosity £, which has often been 
neglected in hydrodynamic simulation of nuclear collisions. The zero bulk viscosity is for a conformal equation of 
state and also a reasonable approximation for the weakly interacting gas of quarks and gluons. For example, the 
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perturbative QCD calculation gives C/s = Q.02a 2 s for 0.06 < a s < 0.3 Q. However, recent lattice QCD results 
show that the bulk viscosity over entropy density ratio C/s rises dramatically up to the order of 1.0 near the critical 
temperature T c [3, G3]. (There are still some subtle issues to determine the bulk viscosity of QCD thr oug h calculating 
the correlations of the energy-momentum tensor on the lattice, see more detailed discussion in Ref. [171].) The s harp 
peak of bulk viscosity at T c has also been observed in the linear sigma model [l|[ and in the real scalar model [l9j |. 
The increasing tendency of C/s has been shown in a massless pion gas [2(| and in the NJL model below T c [HI]. The 
large bulk viscosity near phase transition is related to the non-conformal equation of state p2l . [23[ , and the correlation 
between the bulk viscosity and the conformal anomaly has been investigated in Ref. j24| . 

The sharp rise of the bulk viscosity will lead to the breakdown of the hydrodynamic approximation around the 
critical temperature. The effect of large bulk viscosity on hadronization and freeze-out processes of QGP created at 
heavy ion collisions has been discussed in Refs. [25|, [26|, [2?], [H| . The authors of Ref. [25| pointed out the possibility 
that a sharp rise of bulk viscosity near phase transition induces an instability in the hydrodynamic flow of the plasma, 
and this mode will blow up and tear the system into droplets. Another scenario is pointed out in Ref. (l5l . |27| that 
the large bulk viscosity near phase transition might induce "soft statistical hadronization" , i.e. the expansion of QCD 
matter close to the phase transition is accompanied by the production of many soft partons, which may be manifested 
through both a decrease of the average transverse momentum of the resulting particles and an increase in the total 
particle multiplicity. 

Due to the complexity of QCD in the regime of strong coupling, results on hot quark matter from lattice calculation 
and hydrodynamic simulation are still lack of analytic understanding. In recent years, the anti-de Sitter/conformal 
field theory (AdS/CFT) correspondence has generated enormous interest in using thermal Af = 4 super- Yang-Mills 
theory (SYM) to understand sQGP. The shear viscosity to entropy density ratio rj/s is as small as 1/An in the strongly 
coupled SYM plasma Q. However, a conspicuous shortcoming of this approach is the conformality of SYM: the square 
of the speed of sound c 2 s always equals to 1/3 and the bulk viscosity is always zero at all temperatures in this theory. 
Though C/s at T r is non-zero for a class of black hole solutions resembling the equation of state of QCD, the magnitude 
is less than 0.1 [23], which is too small comparing with lattice QCD results. 

An alternative nonperturbative approach to study QCD phase transition is by using effective models. There is still 
no satisfied dynamic model which can describe deconfinement phase transition successfully. We thus focus on effective 
models which can describe chiral symmetry restoration. QCD with two-flavor massless quarks has a global symmetry 
SU(2)r x SU(2)l, which is isomorphic to 0(4). The chiral symmetry restoration in the case of 2-flavor QCD is of 
second order phase transition, and the universal critical behavior falls in the same universality class as the 0(4) model 
||. It has been argued that the chiral phase transition at finite chemical potential near critical point belongs to the 
Z(2) universality class [J]. This motivates us to study the thermodynamic and transport properties of Z(2) and 0(4) 
models near phase transition in this paper. 

The critical phenomena in the Z{2) and 0(4) models have been well-known, and the singular behavior of the 
static and dynamic properties has been very well studied. However, due to the finite size and time effects for 
the system created in heavy ion collisions, the critical singularity will not show up in the observables. Therefore, 
in this paper, we will not focus on the singularities at the critical point. We will study the thermodynamic and 
transport properties in the Hartree approximation of the Cornwall-Jackiw-Tomboulis (CJT) formalism [301 ] . At finite 
temperature, the naive perturbative expansion in powers of the coupling constant breaks down, and CJT formalism 
provides a convenient resummation method. In the Hartree approximation, we cannot perform the precise critical 
behavior at phase transition, but we can get the qualitative behavior near phase transition. 

It has been found in Ref. [3l| that in the simplest real scalar model with Z{2) symmetry breaking in the vacuum, 
r\/ s behaves the same way as that in systems of water, helium and nitrogen in first-, second-order phase transitions 
and crossover [32] ]. In Ref. [l{|, we have investigated the equation of state and bulk viscosity in the real scalar 
model in the case of 2nd order phase transition, and we have found that the thermodynamic properties and transport 
properties in this simple model at strong coupling are similar to those of the complex QCD system. In this paper, 
we will systematically investigate the thermodynamic properties and bulk viscosity of the Z(2) model and the 0(4) 
model in the cases of first-, second-order phase transitions and crossover in the framework of CJT formalism. 

This paper is organized as follows. In Sec. [UJ we introduce the CJT formalism in the Z(2) and 0(4) model. In 
Sec. IIII1 we introduce the thermodynamic quantities and the bulk viscosity. In Sec. HVi we present our numerical 
results. At the end, we give discussions and summary in Sec. [V] 
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II. Z(2) AND 0(4) MODELS IN THE FRAMEWORK OF CJT FORMALISM 

A. CJT formalism 



At finite temperature T, the temperature introduces a new energy scale which can conspire with the typical 
momentum scale p of a process so that gT/p is no longer of order g, but can be of order 1 [33| . Consequently, all 
terms of order gT/p have to be taken into account which requires the resummation of certain classes of diagrams. The 
CJT formalism, which is equivalent to the ^-functional approach of Luttingcr and Ward [34| and Baym [35| , provides 
such a convenient resummation method. In this paper, we follow the notations used in Ref. [36j ] - 

The CJT formalism can be viewed as a prescription for computing the effective action of a given theory, it generalizes 
the concept of the effective action T[<f>] for the expectation value <p of the one-point function in the presence of external 
sources to that for the effective action r[0,G] for </> and the expectation value G of the two-point function in the 
presence of external sources, with 

Y[4>,G\ =vS] + \ Tr In G" 1 + h^{G^G-l) + T 2 [$,Gt\. (1) 

Here, Tq[4>] is the tree- level action, Gq 1 the inverse tree- level two-point function, and T 2 [(f>, G] the sum of all two- 
particle irreducible (2PI) vacuum diagrams with internal lines given by G. 
The stationary points of this functional, 



smc] 
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= , (2) 

i =l p,G=G 



provide self-consistent equations for the expectation values of the one- and two-point functions (f> and G in the absence 
of external sources, denoted as <p and G, respectively The stationarity conditions Eq. ([2]) for the effective action 
are Dyson-Schwinger equations for the one- and two-point Green's functions of the theory. The Dyson-Schwinger 
equation for the two-point Green's function can be derived and has the form of 

G-^G^+li, (3) 

where 

ST 2 $,G] 



n = -2 



SG 



(4) 

=ip,G=G 



is the self-energy. 

In general, the CJT formalism resums one-particle irreducible diagrams to all orders. As long as I^ contains all 
2PI diagrams, the CJT effective action is exact. However, it is practically impossible to compute all 2PI diagrams, 
and one has to truncate T2 at some order in the number of loops. The advantage of the CJT formalism is that any 
truncation of Y 2 yields a many-body approximation scheme which preserves the symmetries of the tree-level action. 
The solution of Eqs. ^ is thcrmodynamically consistent and conserves the Noether currents. 

The CJT formalism is quite useful for studying theories with spontaneously broken symmetries. In the following, 
we use the CJT formalism studying the Z(2) and 0(4) models with spontaneous symmetry breaking in the vacuum 
and symmetry restoration at finite temperature. 



B. Z(2) model in the CJT formalism 

The critical end point has been studied by using the Ginzburg-Landau effective potential of the order parameter 
field (scalar field) up to the sixth order, e.g. see Ref. [3|. Here we introduce the real scalar theory including the sextet 
interaction which described by the Lagrangian 

£ = \(d^f - \a4? - \h* - \ccf + Hep. (5) 

When H = 0, this theory is invariant under <p — ► —<p and has a Z 2 symmetry. 

Unlike in the Ginzburg-Landau effective potential Q where a, b, c are functions of temperature, here a, b, c are model 
parameters, which determine the vacuum properties. The system at finite temperature will be evaluated in the CJT 
formalism. We will discuss the following four cases: 1) c = 0, b > 0, a > 0, H = 0, the system is always in the symmetric 
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phase. 2) c = 0, b > 0, a < 0, H = 0, the vacuum at T = breaks the symmetry spontaneously, and the symmetry 
is restored at higher T with a second-order phase transition. 3) c = 0, b > 0, a < 0, H ^ 0, the Z(2) symmetry is 
explicitly broken, and the system will experience a crossover at high temperature. 4) c > 0, b < 0, a > 0, H = 0, the 
broken symmetry is restored at high T with a first-order phase transition. 

Assuming translation invariance, we consider effective potential Q, instead of effective action T, these two quantities 
are related via: 



V 
T' 



0, 



where V is the 3-volume of the system. The effective potential in the CJT formalism reads [36 



K 



where 



H<j> 



is the tree-level potential, and G(Gq) is the full(trec-lcvel) propagator: 

G- 1 {K,4>) = -K 2 + m\4>) ,Gv\K,4>) = -K 2 + ml{4>) , 
with the tree-level mass tUq = a + 3b <jr + 5c 4> 4 . In Hartrce approximation, the 2PI potential ^ only includes 



(G) 



(7) 



(8) 
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^G]=[-b+ 1 4i 



K 



K 



The self-consistent one- and two-point Green's functions satisfy 
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This allows us to solve 4> and m through the coupled Dyson-Schwinger equations: 



(10) 



111) 



a<t> + bej) 6 + c<jf + (3b<j) + 10a/) 6 ) / G{K) + 15c0 / G(/\) = H , 



K 



K 



m 2 - mo = 3(6 + 1Oc0 2 ) J G{K) + 15c I J G{K) 



(12) 



C. 0(4) model in the CJT formalism 

The Lagrangian of the O(N) model reads 

£(0) = ia A1 0-5^-^0-0-^ F (0-0) 2 + i?0 1 , (13) 

where (ft = (0i, ...,^>jv) is an O(N) vector. We identify the first component <f>i with the a field and the remaining 
— 1 components as the 7r fields. The last term Htfii breaks the symmetry explicitly and has been introduced in 
order to generate masses for the pions. For H = 0, a > and b > 0, the Lagrangian is invariant under O(N) rotations 
of the fields. For H = 0, a < and b > 0, this symmetry is spontaneously broken down to 0(N — 1), leading 
to N — 1 Goldstonc bosons (the pions), and the field <j> obtains a non- vanishing vacuum expectation value <fi. The 
O(N) symmetry will be restored at finite temperature with a second-order phase transition. By shifting the field as 
4> — > 4> + 4>, the "classical potential" takes the form 



(14) 
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The inverse tree-level propagator which corresponds to the above Lagrangian density is 

where Gq^ , G^ 1 are sigma and pion propagator respectively. 

The CJT effective potential of the O(N) model can be written as the following function of full propagators: 



(15) 
(16) 



n&G^G*) = n o (0) + - J [inG^ik) + G^(k-J)G a (k) - 1} 



+- 



N -1 



[lnG; 1 ^) +G - 1 (fc;^)G 7r (fc)-l] + Q 2 (<£,G CTJ G ff ) 



In the Hartree approximation, the 2PI effective potential takes the form of 
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Jq 


+ 3 4iV 


Jq 



n 2 = (N + i)(N-i) — 



The stationarity condition is written as 

sn[4>, G a , GJ 



+ nN-l)^J Q G a {Q)J L G,{L) 



SQ[(j), Gg, G^} 

SG a {k) 

Sfl[4>, G ai Gtt] 
SG n {k) 



i=<4, G a =G„.G Jr =G 1 , 



=<f>,G tr =G tT ,G„=G 7 , 



= 



o , 



(17) 

(18) 

(19) 
(20) 
(21) 



which determines the expectation values of the one- and two-point functions in the absence of external sources </> and 
G(j , Gjr. Similarly, the Schwinger-Dyson equations for sigma and pion propagators for the effective potential Q are 
given by 



G-\k) = G^{k-4>) + ^{k) = -k 2 + M%, 
G~ l {k) = G^(k; <f) + E.(fc) = -k 2 + Ml . 

By solving Eq. ([2D|) and ([2"Tj) . the corresponding self-energy can be obtained as 



E a (k) 



5G a {k) 



< j > =^ >l G,,=G !T ,G 7 ,=G 7 , 



2 6Q 2 [^,G a ,G v ] 



N-l 5Gn{k) 



(22) 
(23) 



(24) 
(25) 



b=<l>,G <J =G<,,Gn=G* 



In our following numerical calculations, we will take the case of N = 4. 



III. THE EQUATION OF STATE, SOUND VELOCITY, SPECIFIC HEAT, INTERACTION MEASURE 

AND BULK VISCOSITY 



Once we have the effective potential f2, we can derive all thermodynamic properties of the system. The entropy 
density is determined by taking the derivative of effective potential with respect to the temperature, i.e, 

s = -dVL{4>)/dT. (26) 

In the symmetry breaking case, the vacuum effective potential or the vacuum energy density is negative, i.e, 



n v = n(<f>)\ T = < 0. 



(27) 



6 



As the standard treatment in lattice calculation, wc introduce the normalized pressure density px and energy density 
ex as 

p T = -Q, T , £t = -Pr + Ts, (28) 

with 

n T = - n v . (29) 

The equation of state Pt{^t) is arL important input into hydrodynamics. The square of the speed of sound C 2 is 
related to px /ct an d has the form of 

C ' = Te = Td^dT = (30) 

where 

C v = de/dT, (31) 

is the specific heat. At the critical temperature, the entropy density as well as the energy density change most quickly 
with temperature, thus one expect that C 2 should have a minimum at T c . The trace anomaly of the energy-momentum 
tensor T» u 

(32) 

is a dimensionless quantity, which is also called the " interaction measure" . 

The bulk viscosity is related to the correlation function of the trace of the energy- momentum tensor 

1 1 f°° f 

C = - Urn - / dt d 3 re Mt ([OWx), 0£(O)]> . (33) 



9 oj 







Accordingto the result derived from low energy theorem, in the low frequency region, the bulk viscosity takes the 
form of [15| 



9.o ^^t^ + ^I 

= -^—{-16e T + 9TS + TC v + 16\e v \}. (34) 
9 ujq 

with the negative vacuum energy density e v = Cl v = £l(<fr)\T=o, an d the parameter loq = loq(T) is a scale at which the 
perturbation theory becomes valid. From the above formula, we can see that the bulk viscosity is proportional to the 
specific heat C v near phase transition, thus £/s behaves as near T c in this appproximation. 

IV. THERMODYNAMIC PROPERTIES OF Z(2) AND 0(4) MODELS 
In this section, we show our numerical results for thermodynamic properties of the Z(2) and 0(4) models. 

A. Z(2) model without symmetry breaking in the vacuum 

We firstly consider the real scalar model without symmetry breaking in the vacuum, i.e. a > 0. In Fig. [1] (a) — (e), 
we show the ratio of the pressure density over energy density pt/^t, the trace anomaly (ey — ^Pt)/T 4 , the square 
of sound velocity C^, the specific heat C v and the bulk viscosity over entropy density ratio £/s as functions of the 
temperature T for different coupling strength b. The parameters taken for calculations are: 1) a = 1000 McV 2 , b = 0.1, 
2) a = 10000 MeV 2 , b = 10, 3) a = 10000 MeV 2 , b = 30, and 4) a = 10000 MeV 2 , b = 60. 

In the weak coupling case when b = 0.1, it is found that the pressure density over energy density pt/ct, the sound 
velocity square C 2 and the specific heat C v show similar behavior as in the ideal gas. Both pt/^t and C 2 increase 
with the temperature T, and reach the conformal value 1/3 at high temperature. The specific heat C v monotonically 
increases with temperature. There is no much information about the trace anomaly (ey — 3pt)/T 4 for an ideal gas 
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in the literature. It is found that the trace anomaly (ex — 3pr)/T 4 shows a peak at low temperature, then decreases 
monotonically and reaches the conformal value at high temperature. 

The peak of the trace anomaly (ex ~ 3pt)/T 4 at low temperature, which is not related to the phase transition was 
also observed in Ref. [l9| in the real scalar model with 2nd-order phase transition and in Ref. [24[ in the Chiral 
Perturbation Theory for the pion gas. In Ref. (24[, the low-T peak of the trace anomaly was interpreted as the explicit 
conformal breaking, whose contribution comes from massive pions. However, for the real scalar system, there are no 
massive pions, it is not clear for us what is the reason inducing the conformal symmetry breaking at low temperature. 

In the case of strong coupling, we observe that the pressure density over energy density pr/f-T, the sound velocity 
square C 2 increase with temperature and saturate at high temperature. The pressure density over energy density 
Pt/^t saturates at a value smaller than 1/3, the stronger the coupling strength is, the smaller value Pt/^t saturates. 
The sound velocity square C 2 still saturates at 1/3. The specific heat C v increases with temperature. However, it is 
found that the low-T peak of the trace anomaly slowly disappears with the increase of the coupling strength, and at 
high T, the trace anomaly goes to a larger value for stronger coupling strength 6, which indicates that the conformal 
symmetry is broken at high T due to the bare strong coupling. Therefore, the trace anomaly at high temperature can 
indicate the strength of the coupling, in this sense, the trace anomaly is also called the "interaction measure". 

In both weak coupling and strong coupling, the bulk viscosity over entropy density £/s decreases monotonically 
with the increase of the temperature. At high T, £/s reaches its conformal value in the case of weak coupling, and 
reaches a finite value in the case of strong coupling. However, we don't observe the correlation between the trace 
anomaly and the bulk viscosity as shown in Ref. [2J], where it shows that the peak at low temperature in the trace 
anomaly also appears in the bulk viscosity. 

B. Z(2) model with a second-order phase transition 

In Fig. 0(a) — (e), we show the ratio of pressure density over energy density pt/^t, the trace anomaly (er — 3pr)/T 4 , 
the sound velocity square C 2 , the specific heat C v and the bulk viscosity over entropy density (/s as functions of 
the temperature T for different coupling strength b for the real scalar model with a second order phase transi- 
tion. The parameters used for calculations and the corresponding critical temperatures are: l)a = — 100 MeV 2 , b = 
0.1, T c = 88MeV, 2)a = -10000 MeV 2 , b = 10, T c = 82MeV, 3)a = -50000 MeV 2 , b = 30, = 124MeV, 
4) a = -100000 MeV 2 , b = 60, T c = 144 MeV. The main results in this case have been shown in Ref. [l<|. 

In the weak coupling case when b = 0.1, it is found that the pressure density over energy density pt/^t-, the sound 
velocity square C 2 and the specific heat C v show similar behavior as in the ideal gas except near the phase transition 
region. Both pt/^t and C 2 show a downward cusp at T c , and reach the conformal value 1/3 at high temperature. 
The specific heat C v exhibits a small upward cusp at T c . In the weak coupling case, we observe double peak in the 
trace anomaly, one smooth peak shows up at low temperature, another upward cusp appears at T c . 

In the case of strong coupling, at high temperature region T > T c , the behavior of the the pressure density over 
energy density px/er, the sound velocity square C 2 , the specific heat C v , and the trace anomaly show similar behavior 
as those in the symmetric case, pt/^t and C 2 increase with temperature and saturate at high temperature. The 
pressure density over energy density pt/^t saturates at a value smaller than 1/3, the stronger the coupling strength 
is, the smaller saturation value pt/ct will be. The sound velocity square C 2 still saturates at 1/3. The specific heat 
C v increases with temperature. The trace anomaly decreases with temperature and goes to a larger value for stronger 
coupling strength b. 

Near phase transition region T ~ T c , both the pressure density over energy density pt/^t and the sound velocity 
square C 2 show a downward cusp at T c , the specific heat C v and the trace anomaly show an upward cusp at T c . When 
the coupling strength increases, the depth of the downward cusp for the pressure density over energy density pt/ct 
and the sound velocity square C 2 become deeper and deeper, while the height of the upward cusp for the specific heat 
C v and the trace anomaly becomes higher and higher. 

At low temperature region T < T c , both the pressure density over energy density pt/^t and the sound velocity 
square C 2 show a bump, i.e. pt/^t and C 2 firstly increase with T then decrease with T. However, it is found that 
the low-T peak of the trace anomaly slowly disappears with the increase of the coupling strength. 

The bulk viscosity over entropy density £/s decreases with T at low temperature region, then rises up at the critical 
temperature T c , and decreases further in the temperature T > T c . 

C. Z(2) model with explicit symmetry breaking 

In Fig. [3] (a) — (e), we show the ratio of pressure density over energy density pt/ct, the sound velocity square C 2 , 
the trace anomaly (ey — 3pt)/T , the specific heat C v and the bulk viscosity over entropy density (/s as functions of 
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FIG. 1: The ratio of pressure density over energy density pr/er, the trace anomaly (er — 3pr)/T 4 , the sound velocity square 
Cs, the specific heat C v and the bulk viscosity over entropy density ratio (/s as functions of the temperature T for the real 
scalar model in the symmetric phase. 
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FIG. 2: The ratio of pressure density over energy density pr/tr, the trace anomaly (er — 3pr)/T 4 , the sound velocity square , 
the specific heat C v and the bulk viscosity over entropy density ratio £/s as functions of the temperature T with a second-order 
phase transition in the real scalar model. 
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the temperature T for second order phase transition and crossover for the real scalar model. The parameters used for 
calculations are: 1) a = -10000 MeV 2 , b = 10, H = 0, and 2) a = -10000 Me V 2 , b = 10, H = (40MeV) 3 . The only 
difference of the two set of parameters is the value of H. When H = 0, the system experiences a second order phase 
transition, and when H ^ 0, the system experiences a crossover. 

At both low temperature region T <T C and high temperature region T > T c , the behavior of all the thermodynamic 
quantities and bulk viscosity don't show much difference for these two set of parameters. However, near critical 
temperature region T ~ T c , it is observed all the cusp behaviors are washed out, e.g. the downward cusp in the 
pressure density over energy density pt / £t and the sound velocity square C 2 develops into a shallow valley, and the 
upward cusp in the trace anomaly develops into a smooth peak. Especially, the upward cusp for the specific heat C v 
and the bulk viscosity over entropy density f/s vanishes and there is no obvious change near critical temperature. 

D. Z(2) model with a first-order phase transition 

In Fig. 2] (a) — (e) , we show the ratio of pressure density over energy density pt /ct, the sound velocity square C 2 , the 
trace anomaly (ey — 3pt)/T 4 , the specific heat C v and the bulk viscosity over entropy density £/s as functions of the 
temperature T for the real scalar model with a first order phase transition. The parameters used for calculations and 
corresponding critical temperatures are: 1) a = 100 MeV 2 , b = -0.125, c = 0.000025 MeV -2 , T c = 59 MeV, 2) a = 
10000 MeV 2 , b = -1.2, c = 0.000025 MeV" 2 , T c = 155 MeV, 3) a = 40000 MeV 2 , b = -A, c = 0.000069 MeV" 2 , T c = 
204 MeV. 

At high temperature region T > T c , the behavior of the pressure density over energy density pt/ct, the sound 
velocity square C 2 , the specific heat C„, and the trace anomaly show similar behavior as those in the case of symmetric 
phase, pt/^t and C 2 increase with temperature and saturate at the conformal value 1/3 at high temperature. The 
specific heat C v increases with temperature and the trace anomaly decreases with temperature and goes to the 
conformal value 0. 

At low temperature region T < T c , the behavior of the pressure density over energy density Pt/^t, the sound 
velocity square C 2 , the specific heat C v , and the trace anomaly also show similar behavior as those in the case of 
symmetric phase. Both the pressure density over energy density pt/ct and the sound velocity square C 2 monotonically 
increase with the temperature. The low-T peak of the trace anomaly still shows up in the weak coupling case and 
slowly disappears with the increase of the coupling strength. 

Near phase transition region T ~ T c , the behavior of the pressure density over energy density pt/ct, the sound 
velocity square C 2 , the specific heat C„, and the trace anomaly also show similar behavior as those in the case of 
second order phase transition. The only difference is that the width of the cusp becomes very narrow. The downward 
cusp in the pressure density over energy density pt/^t and the sound velocity square C 2 becomes a dip at T c , and the 
upward cusp of the specific heat C v and the trace anomaly develops into a delta function at T c . When the coupling 
strength increases, the value of the pressure density over energy density pt/ct and the sound velocity square C 2 at 
T c become smaller and smaller, while the value of the specific heat C v and the trace anomaly becomes bigger and 
bigger. 

The bulk viscosity over entropy density £/s decreases with T at low temperature region, then sharply rises up at 
the critical temperature T c , and decreases further in the temperature T > T c . The value of £/s becomes larger and 
larger at the critical temperature with the increase of coupling strength. 

E. 0(4) model with spontaneous symmetry breaking in the vacuum 

In Fig. 0a — e, we show the ratio of pressure density over energy density pt/^t, the sound velocity square C 2 , the 
interaction measure (ey — 3j»t)/T 4 , the specific heat C v and the bulk viscosity over entropy density £/s as functions of 
the temperature T for the 0(4) model with a second order phase transition in the chiral limit H = 0. The parameters 
used for calculation arc taken from [37j]: 1) a = -(282.84MeV) 2 , b = 39.53, 2) a = -(424.264MeV) 2 , b = 88.88, and 
3) a = — (565.685MeV) 2 , b = 158.02 which produce the vacuum pion mass = 0, vacuum pion decay constant 
/ w = 90 MeV, and vacuum sigma meson mass as m a = 400MeV m a = 600MeV and m a = 800MeV, respectively. The 
stronger coupling strength of b corresponds to the larger sigma mass in the vacuum. The critical temperature for the 
chiral symmetry restoration in these three cases are T c = 176MeV. 

It is found that all the thermodynamic properties and bulk viscosity show similar behaviors as those in the Z{2) 
model with second order phase transition at strong coupling. 

At high temperature region T > T c , the behavior of the the pressure density over energy density pt/^t, the 
specific heat C„, and the trace anomaly show similar behavior as those in the symmetric case, pt/^t increases with 
temperature and saturate at a value smaller than 1/3, the stronger the coupling strength is, the smaller saturation 
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FIG. 3: The ratio of pressure density over energy density pr/er, the trace anomaly (er — 3pr)/T 4 , the sound velocity square 
Cs , the specific heat C v and the bulk viscosity over entropy density ratio C,/s as functions of the temperature T with a crossover 
and with a second-order phase transition in the real scalar model. 
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FIG. 4: The ratio of pressure density over energy density pr/er, the trace anomaly (er — 3pr)/T 4 , the sound velocity square 
Cs , the specific heat C v and the bulk viscosity over entropy density ratio (/s as functions of the temperature T with a first-order 
phase transition in the real scalar model. 
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value pt/^t will be. The specific heat C v increases with temperature. The trace anomaly decreases with temperature 
and goes to a larger value for stronger coupling strength b. However, it is found that the sound velocity square C 2 
shows different behavior comparing with that in the Z(2) case, C 2 decreases with the temperature in the region 
T > T c in the 0(4) model and saturates at a value smaller than 1/3. 

At low temperature region T < T c , both the pressure density over energy density pt/&t and the sound velocity 
square C 2 show a bump. There is no low-T peak showing up in the trace anomaly, because all these three sets of 
parameters correspond to strong coupling. 

Near phase transition region T ~ T c , both the pressure density over energy density pt/^t and the sound velocity 
square C 2 show a downward-cusp at T c , the specific heat C v and the trace anomaly show an upward-cusp at T c . When 
the coupling strength increases, the corresponding values of the pressure density over energy density Pt/^t and the 
sound velocity square C 2 at T c become smaller and smaller, while the critical values of the specific heat C v and the 
trace anomaly become larger and larger. 

The bulk viscosity over entropy density £/s decreases with T at low temperature region, then rises up at the critical 
temperature T c . The critical value of (/s at T c in the 0(4) model which simulates the real QCD chiral phase transition 
is around 0.02 — 0.06, which is rather small comparing with the lattice QCD results. 

F. 0(4) model with explicit symmetry breaking 

In Fig. [B]a — e, we show the ratio of pressure density over energy density pt/^t, the square of sound velocity O 2 , 
the trace anomaly (er — 3pt)/T 4 , the specific heat C v and the bulk viscosity over entropy density (/s as functions of 
the temperature T for the 0(4) model in the case of crossover when chiral symmetry is explicitly broken by a finite 
value H = (121.6MeV) 3 . The parameters used for calculation arc taken from [37[: 1) a = -(225.41 MeV) 2 , b = 32.12, 
2) a = -(388.34 MeV) 2 , b = 76.17, and 3) a = -(539.27MeV) 2 , b = 145.02, which produce the vacuum pion mass 
to?,- = 139.5 MeV, the vacuum pion decay constant / w = 92.4 MeV, and vacuum sigma meson mass as m a = 400MeV 
m a = 600MeV and m a = 800MeV, respectively. 

At both low temperature region T < T c and high temperature region T > T c , the behavior of all the thermodynamic 
quantities and bulk viscosity over entropy density ratio show similar behavior as the case of second order phase 
transition in the 0(4) model. However, near critical temperature region T ~ T c , it is observed all the cusp behaviors 
are washed out, e.g. the downward cusp in the pressure density over energy density pt/^t and the sound velocity 
square C 2 develops into a shallow valley, and the upward cusp in the specific heat and the trace anomaly develops 
into a smooth peak. The sharp rise of the bulk viscosity over entropy density Q/ s in the case H = vanishes in the 
chiral symmetry explicitly breaking case and there is no obvious change of £/s near the critical temperature. Together 
with the observation in the Z(2) model, we can predict that at RHIC, where there is no real phase transition and the 
system experiences a crossover, the bulk viscosity over entropy density is small, and it will not affect too much on 
hadronization. 



V. DISCUSSIONS AND SUMMARY 
A. Comparing with other results in AdS/CFT, lattice QCD and effective QCD models 

We have investigated the thermodynamic properties and bulk viscosity in the Z(2) and 0(4) models in the Hartree 
approximation of CJT formalism. We now compare our results with the results in AdS/CFT, lattice QCD and effective 
QCD models like the Polyakov-loop Nambu-Jona-Lasinio (PNJL) model. 

The conformal limit has attracted much attention in recent years, since people are trying to understand strongly 
interacting quark-gluon plasma by using AdS/CFT techniques. In conformal field theories including free field theory, 
the pressure density over energy density and the sound velocity square is always 1/3, i.e. pt/^t — <? s = V^j and the 
trace anomaly and the bulk viscosity is always zero, i.e. A = ( = 0. Lattice results show that at asymptotically high 
temperature, the hot quark-gluon system is close to a conformal and free ideal gas. Our results of Z(2) and 0(4) 
models in the weak coupling also reach the conformal limit at high temperature. 

However, lattice results show that near deconfincmcnt phase transition, the hot quark-gluon system deviates far 
away from conformality. Both Pt/ct and c 2 show a minimum around 0.07, which is much smaller than 1/3. For the 
SU(3) pure gluon system[22||, the peak value of the trace anomaly A£^ T reads 3 ~ 4 at T rnax and the corresponding 
"interaction measure" is A^ T /dc = 0.2 ~ 0.25, with the gluon degeneracy factor da = 16. (Note that here 
T max ~ 1.1T C is the temperature corresponding to the sharp peak of A.) For the two-flavor case [23||, the lattice result 
of the peak value of the trace anomaly A^^ 2 reads 8 ~ 11, the corresponding interaction measure at T max is given 



14 




FIG. 5: The ratio of pressure density over energy density pr/er, the trace anomaly (er — 3pr)/T 4 , the sound velocity square 
Cs , the specific heat C v and the bulk viscosity over entropy density ratio f/s as functions of the temperature T with second 
order phase transition in the 0(4) model. 
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FIG. 6: The ratio of pressure density over energy density pr/er, the trace anomaly (er — 3pr)/T 4 , the sound velocity square 
Cs , the specific heat C v and the bulk viscosity over entropy density ratio £/s as functions of the temperature T with explicit 
symmetry breaking in the 0(4) model. 
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as A^jJ^ 2 1 (do + (Iq) = 0.28 ~ 0.4, with quark degeneracy factor dQ = 12. There have been some efforts trying to 
understand trace anomaly in gluodynamics near and above T c in terms of dimension two gluon condensate and an 
effective "fuzzy" bag model j38f . 

Our results in Z(2) and 0(4) models show that at critical temperature T c , the trace anomaly A, the specific heat 
C v as well as bulk viscosity to entropy density ratio C/s show upward cusp at T c , and their peak values increase 
with the increase of coupling strength. The ratio of pressure density over energy density pt/ct and the square of the 
sound velocity show downward cusp at T c , which are similar to the behavior of r\j s found in Ref. [3l| . and the 
cusp values decrease with the increase of coupling strength. These cusp behaviors at phase transition resemble lattice 
QCD results. 

To our surprise, we find that when b = 30, the strongly coulped real scalar system can reproduce all thermodynamic 
and transport properties of hot quark-gluon system near T c . Pt/^t at T c is close to the lattice QCD result 0.07, 
A/d = 0.48 (d = 1 for scalar system) at T c is close to the lattice result of the peak value A^j(y 2 / (da + ^q) — 0.4 at 
T m ax- The bulk viscosity to entropy density ratio (/s at T c is around 0.5 ~ 2.0, which agrees well with the lattice 
result in Ref. [l6j|. (Note, here (/s = 0.5,2 correspond to ui a = 10T, 2.5T, respectively.) More surprisingly, the shear 
viscosity over entropy density ratio r\j s at T c is 0.146, which also beautifully agrees with lattice result 0.1 ~ 0.2 in 
Ref. [Hi . 

Table U shows our results of equation of state and transport properties in the real scalar field theor y o r Z ( 2) model 
at T c for second order phase transition, and corresponding results in lattice QCD calculations [Til. IT5. [22l. l23l| . the 
PNJL model [H,[ll, and black hole duals (2j|. In black hole solutions, it is found that Civ ^ 2(1/3 - C 2 ). In the 
real scalar model near phase transition, there is no any universal relationship between £ and r\. 





v/s 




A/d 


cl 


pr/er 


6 = 30 


0.146 


0.5 ~ 2.0 


0.48 


0.03 


0.07 


LAT G [15, 22] 


0.1 ~ 0.2 


0.5 ~ 2.0 


0.25 




0.07 


LATjv/=2 [15, 23] 




0.25 ~ 1.0 


0.4 


0.05 


0.07 


PNJL [39, 40] 






0.21 


0.08 


0.075 


AdS/CFT [9] 


1/47T 








1/3 


1/3 


TypelBH [29] 


1/47T 


0.06 




0.05 




TypcIIBH [29] 


1/47T 


0.08 




~ 





TABLE I: Thermodynamic and transport properties at T/T c = 1 in Z(2) at strong coupling b = 30, in lattice QCD p~5l. Il6l. I22I. 
|23| |. PNJL model [39. |40||. and black hole dules [2S|. The degeneracy factor d = 1,16,28 for Z(2) model, pure gluon system, 
and 2-flavor quark-gluon system, respectively. Lattice results are taken at Tmax- 



B. Low-T peak of the trace anomaly 

In the Z(2) model, it is found that the trace anomay (er — 3pr)/T 4 shows a peak at low temperature in the case 
of weak coupling, and this low-T peak disappears in the strong coupling case. 

The peak of the trace anomaly (et — ^Pt)/T at low temperature, which is not related to the phase transition, 
was also observed in Ref. [24| in the Chiral Perturbation Theory for the pion gas. In Ref. [24], the low-T peak of 
the trace anomaly was interpreted as the explicit conformal breaking, whose contribution comes from massive pions. 
However, for the real scalar system, there are no massive pions, it is not clear for us what is the reason inducing the 
conformal symmetry breaking at low temperature. 

We don't observe the correlation between the low-T peak of the trace anomaly and the bulk viscosity. It might be 
due to the method we have used to calculate the bulk viscosity. It is worthy of checking whether the correlation really 
exists by using Kubo formula to calculate the bulk viscosity. 

C. C/s at RHIC 

Recent lattice QCD results show that the bulk viscosity over entropy density ratio C/s rises dramatically up to the 
order of 1.0 near the critical temperature T c [lj| [l6[. The sharp rise of the bulk viscosity will lead to the breakdown 
of the hydrodynamic approximation around the critical temperature, and will affect the hadronization and freeze-out 
processes of QGP created at heavy ion collisions. The authors of Ref. (2f| pointed out the possibility that a sharp 
rise of bulk viscosity near phase transition induces an instability in the hydrodynamic flow of the plasma, and this 
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mode will blow up and tear the system into droplets. Another scenario is pointed out in Ref. [15|, [2jj that the large 
bulk viscosity near phase transition might induce "soft statistical hadronization" . 

However, if the strongly coupled QGP created at RHIC experiences a crossover and not a real phase transition, 
from our results in the simple Z(2) and 0(4) models, the sharp rise of the bulk viscosity over entropy density will 
be washed out. In this case, one does not need to worry about the breakdown of the hydrodynamic approximation 
around the critical temperature. 

D. Using C/s to locate the CEP 

For real QCD with two quarks of small mass, it is expected that there exists a critical end point (CEP) in the T — fj, 
QCD phase diagram. At small baryon chemical potential /x, the chiral phase transition is a smooth crossover at finite 
temperature. At finite baryon chemical potential, the chiral phase transition is of first order. The precise location 
of the CEP is still unknown. In the future plan, RHIC is going to lower the energy and trying to locate the CEP. 
Recently, the authors of Ref. @ suggested using the shear viscosity over entropy density ratio 77/s to locate the CEP. 

From Rcfs.[3l], [HI, we know that 77/s shows a shallow valley in the case of crossover and a jump at T c for first- 
order phase transition. But it is hard to distinguish whether the system experiences a crossover or first-order phase 
transition just from the value of 77/s extracted from the elliptic flow u 2 . 

From our results in Z{2) and 0(4) models, it is found that the ratio of Q/s shows a very sharp peak at T c , and 
there is no obvious change of (/s for crossover. As pointed out in Ref. [25[ that a sharp rise of bulk viscosity near 
phase transition induces an instability in the hydrodynamic flow of the plasma, and this mode will blow up and tear 
the system into droplets. Therefore, one can distinguish whether the system experiences a first order phase transition 
or a crossover from observables at RHIC experiments. 

Therefore, Q/s is a better quantity than 77/s to locate the CEP. (It is noticed that because the QGP created at 
RHIC is a finite system, we don't discuss the singularity of rj/s,Q/s at CEP. ) 

E. Limitations of our results 

At the end, we hope to point out some limitations of our results. 

Firstly, the results of thermodynamic properties in this paper are based on Hartrcc approximation in the CJT 
formalism. As wc know that mean-field approximation cannot describe critical phenomena very well. For 2nd-order 
phase transition in Z(2) model, the specific heat C v should diverge at the critical point, and behave as t~ a near the 
critical point, with t = (T — T c )/T c and a = 0.11. However, in Hartree approximation of CJT formalism, though we 
observe the weak divergence of C v at T c in the case of strong coupling, we can only see a weak upward cusp of C v at 
T c in the case of weak coupling. As pointed out in the introduction, the QGP created at RHIC is a finite system, the 
singularity of 77/s, C/s at critical point will not show up in the observables. If we are only interested in the qualitative 
properties near phase transition, the Hartree approximation can give the dominant contributions. 

Secondly, the results of bulk viscosity in this paper are based on Eq. (|54"|) . The limitation of Eq. (|34p has been 
analyzed in Refs. [4l| and [ijj. From Eq. we see that the bulk viscosity is dominated by C v at T c . If C v diverges 
at T c , the bulk viscosity should also be divergent at the critical point and behave as t~ a . However, the detailed 
analysis in the Ising model in Ref. [12] shows a very different divergent behavior £ ~ t~ zu+a , with z ~ 3 the dynamic 
critical exponent and v ~ 0.630 the critical exponent in the Ising system. 

Thirdly, we should keep in mind that our results at strong coupling are from effective theory. From renormalization 
analysis, the scalar theory will hit a Landau pole when j3b > b with the j3 function (3b = 9b 2 /l6ir 2 . The results for 
large b in the scalar theory are not guaranteed to be valid in the CJT formalism. 

More careful study on thermodynamic properties in this model beyond mean-field approximation is needed, and 
the bulk viscosity with full spectral function of the pressure-pressure correlator is in progress. 

F. Summary 

In summary, in the Hartree approximation of CJT formalism, we have investigated the thermodynamic properties 
and transport properties of the Z(2) model and 0(4) model, and compared these properties in the cases of first and 
second order phase transitions, in the case of crossover and the case of symmetric phase. 

We have seen that at phase transition, the system either in weak coupling or strong coupling shows some common 
properties. 1) The pressure density over energy density ratio pr/c-T, the square of the speed of sound C 2 as well as 
77/s exhibit downward cusp behavior at T c . 2) The trace anomaly A, the specific heat C v as well as £/s show upward 
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cusp behavior at T c . The cusp behavior is related to the biggest change rate of entropy density at T c [43[. 3) The cusp 
behavior in the first order phase transition is sharper and narrower than that in the second order phase transition. 4) 
In the case of crossover, the cusp behavior is washed out. 

Therefore, if the strongly coupled QGP created at RHIC experiences a crossover and not a real phase transition, 
one does not need to worry about the breakdown of the hydrodynamic approximation around the critical temperature. 
Because the behavior of C/s is so different in the case of first order phase transition and the crossover, we also suggest 
that (/s is a better quantity than r\j s to locate the CEP. 
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